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ABSTRACT 

We investigate the consequences of applying different star formation laws in the galaxy for- 
mation model GALFORM. Three broad star formation laws are implemented: the empirical 
relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of 
Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against 
observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We 
start from published models, and investigate which observables are sensitive to a change in the 
star formation law, without altering any other model parameters. We show that changing the 
star formation law (i) does not significantly affect either the star formation history of the uni- 
verse or the galaxy luminosity functions in the optical and near-IR, due to an effective balance 
between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents 
of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that 
a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. 
We show that this plane can be used to discriminate between the star formation laws. 
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1 INTRODUCTION 

A proper understanding of how galaxies form and evolve must in- 
clude a description of star formation and the physics that regulate 
this phenomenon. The lack of a theoretical description of star for- 
mation (SF) has forced g alaxy formation model l ers to adopt simple 
parametric recip e s (e. g. White & Fre nklll99ll; lLacev et al. 1993; 
Cole et al. 1994; Kauffmann & Char lotl 1 19981: ICole et alj l2OO0l: 
Springel et al.l 1200 ll; jCoral 120061 : iMonaco et al.|[2007l ; lLagos et"ail 
20081; ICattaneo et alj2008l) . 

Two forms commonly adopted for the SF law relate the star 
formation rate (SFR) per unit are a, Esfr, to the gas surface den- 
sity either as: (i) Esf r oc E^ s jSchmidt| [l959). or (ii) Esfr oc 
Egas/Tdyn (Shu 1973), where Td yn is a dynamical timescale. Ob- 
servational constraints on the form of the SF law remained scarce 
for many years. Early observations of individual galaxies gave con- 
flicting results for the power-law index N in relation (i) above , re- 
turning values o ver the range N « 1 - 3 (see lSanduleakll 19691 and 
lHartwicklll97 ll for extreme examples). Nor was it clear, when ap- 
plying form (ii), whether a local or global dynamical timescale Td yn 
is more relevant (e.g. the loca l free-fall time of the gas or the or- 
bital time around a galaxy, see lMadordfl977l and lSolomon & Sage! 
1 1988L respectively). The situation improved with the observational 
sample constructed by Kennicutt (1998, hereafter K98), which he 
used to examine correlations between the global SFR and gas sur- 
face densities for a wide range of nearby galaxies, from normal 
spirals to starbursts. K98 found that the global relation was re- 
markably well fit by (i), with A^ = 1.40 ± 0.15, but also nearly 



equally well fit by (ii), when taking Td yn equal to the disk orbital 
time. However, observational studies of the resolved radial profiles 
of SFR and gas in individual galaxies showed a breakdown of the 
power-law relation at low gas surface densities, with a sharp cut- 
off in Esfr below a criti cal gas surface density fcenmcutdll989l: 
iMartin & KennicutfeOQll) . This feature was interpreted using sim- 
ple gravitational arguments as marking the onset of dynamical sta- 
bility in the gas layer, thereby allowing the fragmentation of the gas 
into st ar-forming clo uds at gas surface densities above a threshold 
value jToomre|[l964l) . 

Over the past ten years, enormous advances in the charac- 
terisation of the SF law have been made possible through the 
availability of new, high quali ty, spatially res olved observations 
of HI (e.g rwalteretalj|2008h and CO (e.g. iHelfer etai] 120031: 
iLerov et al]|2009h . and of UV and IR SFR tracers in samples of 
nearby galaxies. These data have allowed the accurate estimation 
of the molecular and atomic hydrogen content of galaxies over a 
wide range of morphologies and gas fractions. In a ddition, more re- 
liable estimates of unobscured SF in the UV ( e.g.lGil de Paz et ail 
l2007h and of dust-obscured SF in the IR (e.g. lCalzetti et alj|2007l) 
have allowed better determinations of the SFR, both globally and 
within individual SF regions. There is now support for a correla- 
tion between Esfr and the molecular gas surface density, S mo i, 
of the form Esfr oc E mo i dWong & Blitzll2002l : iKennicutt et al.l 
l2007l : lBigiel" et al.ll2008h . This relation is much stronger than that 
with either the total cold gas or atomic hydrogen surface densi- 
ties. At high redshifts, th ere are indications that the same form , 
Esfr oc E mo i could hold (Bou che et al.l2007l : lGenzel et al.l201oh . 
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The correlation of Esfr with E mo i seems physically reasonable 
since stars are observed to form in dense molecular gas clouds (see 
ISolomon & Vanden Boutll2005l for a review). 

The thresho ld in the SF law suggested by the observations of 
iKennicuttl (1989) and iMartin & K ennicutt < l200lh has since been 
studied in more detail. SF activity has been observed in galax- 
ies even at low gas surface densities (S ga s ^ IOMq pc~ 2 ; 
Hever et alj|2004 iFumagalli & Gavaz zi 2008; Bigiel et all [20081 ; 
Rovchowdhurv et alj 1 20091 ; IWvder et all 120091) . but following a 



steeper relation in E gas compared to the SF law found at higher 
gas surface d ensities, suggesting a differe nt re gime of SF activity . 
Furthermore, iBlitz & Rosolowskvl d2006l) and iLerov etail d2008l) 
found that the ratio of molecular-to-atomic hydrogen surface den- 
sities, R mo \, is well fitted by a power-law function of the hydro- 
static pressure in galactic disks, i? mo i oc P" xt , with the same rela- 
tion holding from dwarf to normal spiral galaxies. An empirical SF 
law emerges from these studies in which the relationship between 
Rmoi an d P cx t naturally divides the Esfr — E gas plane into two 
regimes of SF activity when a linear relation such as Esfr oc E mo i 
is assumed, according to whether i? mo i > 1 or < 1. In this pic- 
ture, it is the transition from atomic- to molecular-dominated gas 
which causes the change of slope in the Esfr — E gas relation, 
rather than the stability of the disk. Note that the empirical SF 
laws discussed above attempt to describe the SF activity in galaxies 
on scales larger than ~ 500 pc, averaging over representative vol- 
umes of the ISM including both star-forming clouds and intercloud 
gas (e .g. lOnodera et alll201ol ; lFumagalli et al.ll2010l ; ISchruba et all 

limb- 

Even though the quality of the observational data has im- 
proved the characterization of the SF law, there is still un- 
certainty regarding the main physical mechanism that governs 
the SFR in galaxies. Possibilities currently under discussion in- 
clude the formation of molecules on dust grains in turbulence- 
regul a ted giant molecular cl ouds (GMCs; iKrumholz & McKeg 
120051 ; IKrumholz et alll2009bh and co l lisions of GMC s in shear- 
ing disks dTaid l200Cl: ISchavd |2004 Isilk & Normanl 120091 ; see 
|McKee_&p striker 2007 for a comprehensive review of the the- 
ory of SF activity). Direct compariso ns of empirical a nd theoret- 
ical SF models with observations by ILerov et al I d2008h revealed 
that th resholds of large scale stability, as proposed by IKennicuttl 
d 1998)1 or a simple dependence on the orbital or free-fall timescales 
as in mini (feoOO), do not offer a good explanation for the observed 

Esfr oc E mo i relation. On the other hand, for a suitable choice of 

I 1 J — — I 

parameters, the Krumholz et al. (2009b) SF law is able to reproduce 

the observational results described above fairly well. The nature of 
the physical processes underlying the global SF law could be iden- 
tified by resolving individual GMCs and SF regions in the Milky 
Way and other galaxies. Ongoing observations with Herschel and 
future telescopes such as ALMA (Atacama Large Millimeter Ar- 
ray) and the SKA (Square Kilometre Array), will help us to better 
characterise the properties of individual SF regions, and to test the 
applicability of locally-derived SF laws at higher redshifts. 

In this paper, we implement new recipes for SF in the 
GALFORM semi-analytical model of gal axy format i on in a A 
Cold Dark Matter (ACDM) cosmology dCole et all |200C|). We 
focu s on the empirical laws proposed by KennicutJ ~ 19981) 
and Blitz & Rosolowskv (2006) and the theoretical model of 
IKrumholz et all d2009bl) . Our aim is to identify the galaxy observ- 
able properties which can distinguish between different SF pre- 
scriptions, and also to identify behaviour which is common to dif- 
ferent SF laws. To do this, we insert the new SF prescriptions into 
existing GALFORM models (Baugh et al. 2005; Bower et alj|2006t) 



without changing any othe r model paramet ers to achieve a new 
'best' fit to observations (see lCole et alj20o"fj| and lBower et alj2010l 
for discussions of fitting model parameters). We first study predic- 
tions for the cosmic SFR evolution and then for the evolution of the 
cold gas content of galaxies and their level of SF activity, which 
are very sensitive to this choice. An important result obtained in 
this work is that the study of the SFR in galaxies of different stel- 
lar masses offers a useful constraint o n the choice of S F law. New 
radio telescopes such as MeerKAT (Boo th et al"]|2009l) . ASKAP 
djohnston etai]|2008l) and SKA dSchilizzi et alj|2008l) will mea- 
sure the abundance of HI in the Universe and its evolution. ALMA 
will soon begin to probe both the molecular gas contents and dust- 
obscured SFRs in high-redshift galaxies. The predictions made in 
this paper can be used in conjunction with these future observations 
to constrain the physical processes that trigger SF in galaxies. 

This paper is organized as follows. In Section[2] we present the 
GALFORM model and the two variants we use to study the impact 
of applying different SF laws, pointing out the main differences and 
similarities between them. We also describe the SF laws tested. In 
Section [3] we study the impact of changing the SF law on the evo- 
lution of the cosmic SFR density. Section[4]examines the evolution 
of the cold gas content of galaxies, and we discuss the physics that 
shapes its evolution. In Section[5] we present a new means to distin- 
guish between different SF laws, namely the shape of the sequences 
in the SFR-stellar mass plane, and compare with available obser- 
vations. Finally, in Section[6] we discuss our results in the context 
of future surveys and observations. The appendices cover the de- 
tails of the numerical scheme used to implement the new SF laws, 
some details of the SF laws used and present some predictions for 
observables which are insensitive to the choice of SF law. 



2 MODELLING THE STAR FORMATION ACTIVITY IN 
GALAXIES 

We first give an overview of the GALFORM galaxy formation 
model (§ 2.1), before describing the star formation recipes used in 
GALFORM (§ 2.2). In § 2.3 we outline the new star formation laws 
(see Appendix |A] for a description of the numerical implementa- 
tion of these laws and Appendix |B]f or some more details of the SF 
laws). 



2.1 The galaxy formation model 

We study the impact on the galaxy population of applying different 
star format ion laws in the GALFORM semi-analytic model of galaxy 



formation dCole et alll200ol;lBenspn et alj2003l;lB augh et al j2005l; 
iBower et all 120061 ; lLacev et all |2008|; iBenson & Bowed l2010h . 
GALFORM models the main physical processes that shape the for- 
mation and evolution of galaxies: (i) the collapse and merging of 
dark matter (DM) halos, (ii) the shock-heating and radiative cooling 
of gas inside DM halos, leading to the formation of galactic disks, 
(iii) quiescent star formation in galactic disks, (iv) feedback from 
supernovae (SNe), from AGN heating and from photo-ionization 
of the IGM, (v) chemical enrichment of stars and gas, and (vi) 
galaxy mergers driven by dynamical friction within common DM 
halos, which can lead to the formation of spheroids an d trigger 
bursts of SF (for reviews see Baugh 2006; Benson 2010). Galaxy 
luminosities are computed from the predicted star formation and 
chemical enrichment histories using a stellar population synthesis 
model. Dust extinction is calculated self-consistently from the gas 
and metal contents of each galaxy and the predicted scale lengths 
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of t he disk and bulge comp onents using a radiative transfer model 
(see lCole et a"i]|2000l and lLacev et"aT1l201 ll) . 

We consider the models published by Baugh et al. (2005; 
hereafter Bau05) and Bower et al. (2006, hereafter Bow06). 
The main successes of the Bow06 model are the reproduction 
of the observed break in the galaxy luminosity function (LF) 
at z — in the bj- and A'-bands, the observational inferred 
evolution of the stellar mass functions up to z « 5, the bi- 
modality in the colour-magnit ude relation and the abun dance 
and properties of red galaxies dGonzalez-Perez et al, I l2009h . The 
Bau05 model likewise matches bj- and if-band LFs at z = 
0, but also the number counts and r edshift distribution o f sub- 
milli metre selected galaxies (see also I Gonzalez et alj|201lh . disk 
sizes dAlmeida et alj|2007 ; Gonzalez et alj 200S ), the evolution of 
LvQ-emitters dLe Delliou et al .l2006r . lOrsi et"ai] 2008). counts, red- 
shif t distributions and LFs o f galaxies selected in the mid- and far- 
IR dLacev et alj|2008l. l201Ch and the Lyman-break galaxy LF at 
redshifts z = 3 - 10 dLacev et alj201lh . 

We refer the reader to the original papers for full de- 
tails of the Bau05 and Bow06 models. Discussio ns of the dif- 
ferenc e s between the model s can be found in lAlmeida et al.1 
d2007h, IGonzalez et alj d2009l) . I Gonzalez-Perez et al.1 d2009h and 
lLacev et alj d201 II) . Briefly, the main differences between the mod- 
els are: 

• The SF timescale in galactic disks. The Bow06 model adopts 
a parametrization for the SF timescale which depends on the global 
dynamical timescale of the disk, whereas the Bau05 model adopts 
a timescale which depends only on the disk circular velocity (see 
§2.1). Consequently, high redshift galaxies in the Bau05 model 
have longer SF timescales than those in the Bow06 model and are 
therefore more gas-rich. This leads to more SF activity in bursts in 
the Bau05 model compared to the Bow06 model. 

• The stellar initial mass functi on (IMF) i s assum ed to be uni- 
versal in the Bow06 model, with the lKennicuttl J 1983b IMF adopted. 
The Bau05 model assumes two IMFs, a top-heavy IMF during star- 
bursts (SB) and a Kennicutt IMF during quiescent SF in disks. 
Defining the IMF slope x through dN/dlnm oc m~ x , the Ken- 
nicutt IMF has x — 0.4 below m = lMg and x — 1.5 at higher 
masses, while the top-heavy IMF has x = 0. 

• The gas expelled from galaxies by SN feedback is assumed to 
be reincorporated into the hot halo after a few dynamical times in 
the Bow06 model. In the Bau05 model, this happens only after the 
DM halo doubles its mass. This time, called the halo lifetime, is 
usually much longer than the dynamical time of the halo. 

• In order to reproduce the bright-end of the LF, the Bau05 
model invokes SNe driven superwinds with a mass ejection rate 
proportional to the SFR (note that the expelled material is not rein- 
corporated into any subsequent halo in the DM merger tree). The 
Bow06 model includes AGN feedback to switch off the cooling 
flow in haloes where the central black hole (BH) has an Edding- 
ton luminosity exceeding a multiple of the cooling luminosity (the 
multiple being a model parameter). 

• T he Bow06 model includes disk instabilities (e.g. IMo et al.1 
Il998h as an extra mechanism to trigge r bursts of SF and to make 
spheroids, and hence to build BH mass iFanidakis et alj2010h . 

We use halo merger trees extracted from the Millennium cos- 
mological N-body simulation dSpringell2005h in the Bow06 model 
and Monte-Carlo trees in the Bau05 model, as in the original mod- 
els. We use the Monte-Carlo algorithm o f Parki nson etaD d2008h 
to generat e DM merger trees for the B au05 model from that of 
ICole et al.1 d"2000h . IGonzalez" et al.1 d201ll) tested that changing the 



tree-generation scheme does not significantly change the model 
predictions. The Bau05 model adopts a ACDM cosmology with 
a present-day matter density parameter, f2 m = 0.3, a cosmolog- 
ical constant, Qa = 0.7, a baryon density, f^baryons — 0.04, a 
Hubble constant h — 0.7, and a power spectrum normalization of 
erg = 0.93. The Bow06 model is set in the Millennium simulation 
which has the following cosmology: Q m — Qdm + ^baryons = 
0.25 (giving a baryon fraction of 0.18), fi A = 0.75, a a = 0.9 
and h — 0.73. The resolution of the Af-body simulation is fixed at 
a halo mass of 1.72 x lO lo /i~ 1 M . In the Bau05 model, Monte- 
Carlo merger histories are generated for a grid of halo masses. The 
grid is laid down at selected redshifts, sampling haloes in a fixed 
dynamic range of mass which scales as (1 + z)~ A to approximately 
track the break in the mass function. The minimum halo mass for 
merger trees at z = is 10 9 hT 1 Mq. The Monte-Carlo trees there- 
fore include lower mass haloes than in the N-body case. 

2.2 The original star formation laws in GALFORM 

The original GALFORM models used a parametric SF law to com- 
pute quiescent SFR at each timestep dCole et all2 000) 
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(1) 



where M co id is the total cold gas mass in the galactic disk and 
is the timescale for the SF activity which may depend on galaxy 
parameters but not on the ga s mass. This cho ice was motivated by 
the observations of K98 (see Bell et al. 2003). The SF timescale in 
Eq.[T]was parametrized bv lCole et alj (120000 as, 

Tdisk 
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(2) 



where Vo = 200kms _1 , Vdisk is the disk circular velocity at 
the half mass radius, Tdisk = fdisk/Vdisk is the disk dynamical 
timescale, and and a* are free parameters, with values origi- 
nally chosen to repro duce the gas-to-l uminosity ratio observed in 
local spiral galaxies dCole et al]|2000h . Bow06 adopted values of 
e* = 0.0029 and a* = —1.5, mainly to reproduce the galaxy LF, 
without particul ar reference to the cold gas to luminosity ratio (see 
iKim et alfeoioh . 

In the Bau05 model a different SF timescale was adopted, in 
order to reduce the amount of quiescent SF activity at high redshift, 
thereby allowing more SF in bursts. The recipe used in the Bau05 
model is 

T*=7tt(Vdtak/V r ) a *, (3) 

where Vo = 200 km s" 1 , to and a are free parameters set to 
ro = 8 Gyr and a* = —3, in order to match the local gas mass to 
luminosity ratio (see lPower et alj2 010). 

During SBs caused by galaxy mergers and, in the case of the 
Bow06 model, also by disk instabilities, the available reservoir of 
cold gas is assumed to be consumed in the SB event with a finite 
duration. The SF timescale in bursts is taken to be /dynTbui ge , with 
a floor value r.burst.mm, with Tb u i gc being the bulge dynamical 
time. Bau05 adopted /d yn = 50 and T*burst,min = 0.2 Gyr while 
Bow06 used / dyn = 2 and T„ b urst,mm = 0.005 Gyr. 

The need for two free parameters in Eqs. [2] and [3] reflects a 
lack of understanding of the physics of SF. New, high resolution, 
spatially resolved data reveal that a SF law of the form of Eq. Q] 
with t* indep endent of the gas density results in a poor fit to the 
observations dLerov et alj|2008h . It is necessary to revisit the SF 
law used in GALFORM and to study more consistent ways to char- 
acterise the SF in galaxies. 
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2.3 The new star formation laws 

Here we summarise the three new forms of SF law that we imple- 
ment for q uiescent SF in GA LFORM. These are the empirical rela - 
tions of (i) lKennicuttl dl998h an d (ii)lBlitz & Rosolowskvl bOOfjl) . 
and (iii) the theoretical model of iKrumholz et ail d2009bl) . The ob- 
servational situation for starbursts is less clear than it is for disks, so 
we retain the original GALFORM prescription for the star formation 
timescale in SB (see above). 



2.3.1 The Kennicutt-Schmidt law 

Following the pioneering work of lSchmidtl dl959h . observational 
studies have found that the surface density of SF (Esfr) correlates 
with the projected gas density (E gas ). K98 fitted the relation 

Ssfr = AE gaji , (4) 

where N = 1.4 ± 0.15 and A = 0.1470 when E gas and Ssfr 
are measured in A1qpc~ 2 and Mq pc _2 Gyr _1 , respectively. This 
relation holds over 5 orders of magnitude in SFR and gas sur- 
face density, but show a break to a steeper relation in the outer 
regions of spirals and in dwarf galaxies (K89, K98, Bigi el et all 
120081 : iLerov etai]|2008h . 

The Kennicutt-Schmidt SF law combines the power-law de- 
pendence of the SFR on E gas at high gas densities with a cut-off in 
SF below a critical gas surfac e density, S cr it, as observe d at low gas 
surface densities by K89 and Martin & Kennicutt (2001). The E cr it 
threshold is motivated by the lToomrd q 19641) stability criterion. In 
the case of a thin isothermal gas disk with a flat rotation curve, the 
critical surface density for gravitational instability of axisymmetric 
perturbations is given by 



Ecrit 



V2 



■it 7rG 



V 

R' 



(5) 



where Q cr it is a dimensionless constant ~ 1 and <r g is the ve- 
locity dispersion of the gas (see Appendix IB II for a derivation of 
Eq.[5). We adopt cr B = lOkms -1 consistent with the observations 
of lLerovetai](l2008l) . K98 found Q crit « 2.5 (after scaling by our 
choice of <r g , which is larger than the original value adopted by 
K98). By adopting these values, we calculate E cr it, an d suppress 
SF at radii in which E gas < E cr it . 

We will consider two forms of the Kennicutt-Schmidt law, 
with and without SF suppression below E cr it (Eq. |5). We denote 
these SF laws as KS. thresh (with suppression) and KS (without 
suppression). 



2.3.2 The Blitz & Rosolowsky model: molecular gas fraction 
determined by pressure 

Blitz & Rosolowsky (2006, hereafter BR06) based their SF 
law on two observationally motivated considerations: (i) Ob- 
servations in the IR and at millimetre wavelengths suggest 
that stars form in dense gas environments, namely GMCs (see 



1 We have rescaled the value of A to account for the difference between the 
Salpeter IMF assumed by K98 and the Kennicutt IMF used in our model for 
quiescent SF. The SFRs used in K98 were inferred from Ha luminosities, 
and for a fixed Ha luminosity, the SFR on assuming a Kennicutt IMF is 
0.91 times that inferred on adopting a Salpeter IMF. Note also that we al- 
ways use gas surface densities including helium, assuming a hydrogen mass 
fraction of 74%. 



ISolomon & Vanden Boutl2005l for a review). BR06 assume that the 
SFR is set by the surface density of molecular gajf] E mo i, with a 
proportionality factor (given as an inverse timescale) vsf'- 



Esfr = ^SF E n 



(6) 



(ii) BR06 propose that the ratio of molecular to atomic hydrogen 
gas, Rm i, is given by a power law in the internal hydrostatic pres- 
sure in a galactic disk, P ex t , 



Rti 



E(H 2 ) 
E(ffl) 



Pext 



(7) 



BR06 found that the observed molecular-to-atomic ratios 
in their galaxy sample were well fitted using values of 
loK(P /k B [cm- 3 K 1) = 4.54 ± 0.07 and a = 0.92 ± 0.07. 
ILerov et Id] ( 120081) found similar values using a somewhat larger 
sample. The hydrostatic pressure in disk galaxies at the mid-plane 
is calculated following Elmegreen (1993; see Appendix |B2] for de- 
tails). 

The original SF law of Eq. [6] can hence be rewritten in terms 
of the total gas surface density and the molecular-to-total hydrogen 
ratio, / mo i = E mo i/E gas = fl mo i/(-Rmoi + 1), as 



= ^SF f, 



mol <^ ga s - 



(8) 



We consider two case s for ^sf; (i) a co nstant value ^sf = Vsf = 
0.525 ± 0.25 Gyr -1 jLerov et alj200#l which gives a linear de- 
pendence of Esfr on E gas at high surface densities; (ii) a surface 
density dependent vsf given by 



o 

^SF 



1 + 



Sgas 

"e7 



(9) 



where Eo = 200 Mq pc -2 and q = 0.4 are chosen to recover the 
K98 law at high gas surface densities, and the steepe ning seen in 
the S sfr - E gas relation when starbursts are included ( Bigi el et al.1 
2008). The surface density S is similar to the typical surface den- 
sities of individual GMCs in spiral galaxies, so the transition to the 
steeper dependence could be interpreted as happening when these 
clouds start to overlap. We refer to the SF law with a constant ^sf as 
BR and to the second form, where z^sf = usf (E gas ), as BR.nonlin. 



2.3.3 The Krumholz, McKee & Tumlinson model: 
turbulence-regulated star formation activity 

Krumholz, McKee & Tumlinson (2009; hereafter KMT09) calcu- 
late ^sf and / mo i of Eq[8]for a spherical cloud with SF regulated 
by supersonic turbulence. 

KMT09 assume / mo i is determined by the balance between 
the photo-dissociation of H2 molecules by the interstellar far-UV 
radiati on and the formation of molecules on the surface of dust 
grains. IKrumholz et alj d2009al) calculated / mo t theoretically for 
spherical cloud, and showed that it is approximately a function of 
the gas surface density of the atomic-molecul ar complex (or GMC) 
Ecomp and of the gas metallicity Z (see also lMcKee & Krumholz! 
l20ld) . We use eqn.(2) from KMT09 for / mo t (Ecomp, Z). KMT09 
assume that the surface density of the GMC is related to that of the 



2 Note that we include the associated helium in the molecular gas mass, in 
th e same way as for t he total gas mass. 

3 ILerov et all 120081) derived their SFRs from a calibration using far-UV 
luminosities, assuming a lKroupal Jjjoil) IMF. For the Kennicutt IMF, we 
infer SFRs 1 .02 times larger. For simplicity, we do not apply this correction 
factor, since it is so close to 1 . 
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ISM on larger scales by E CO m P = cE gaa , where the clumping factor 
c 1 is a free parameter, which is not predicted by the model. It is 
implicitly assumed that the fraction of the ISM in GMCs is equal to 
the molecular fraction within a single GMC. Given that the KMT09 
formula predicts that / mo j — > in very metal poor environments, 
we assume that a minimum = 10~ 4 characterises pristine gas 
at very high redshift (see Appendix IB 3 . 1 1 for more details). 

KMT09 define ^sf as the inverse of the timescale required 
to convert all of the gas in a cloud into stars, and obtain it from a 
theor etical model of turbulent fragmentation iKrumholz & McKeel 
l2005t) . The parameter vsv in a cloud depends on the cloud surface 
density, E c ;. KMT09 assume that E c ; is a constant (Eo) in normal 
spiral galaxies, but increases in higher density environments, so that 
E c ; = max[Eo, E gas ] (see Appendix lB3.2l for more details). Com- 
bining all of these ingredients, KMT09 obtain: 

/y, \ -0.33 

^ SF (E gas ) = „<> F xf^pJ for E g as<E (10) 

/y \ 0.33 

= ^SFxf-grpJ for E gas >E (11) 

with v$ F = 0.38 Gyr" 1 and E = 85 M Q pc" 2 . The molecular 
fraction / mo i(cE gas , Z) contains the dependence on the clumping 
factor c. The KMT09 SF model thus predicts three regimes of SF 
(qualitatively similar to the BR.nonlin model): (i) regions with low 
gas surface densities and hence low molecular gas fractions, which 
produce steeper SF relations; (ii) an intermediate density regime 
characterised by constant GMC surface densities; (iii) a high den- 
sity regime where there is an increase in the GMC density. This 
model is able to reproduce reasonably well the observed trend of 
Esfr with E gas in nearby galaxies on adopting c « 5, roughly 
consistent with observed overdensities of molec ular gas complexes 
dRosolowskv & Blitzl20"0^ ; ISchuster et alj2007l) . Hereafter we will 
refer to this SF law as KMT. 



3 THE EVOLUTION OF THE SFR DENSITY 

Our aim is to investigate the impact on galaxy properties of chang- 
ing the SF law over a wide range of redshifts, and to compare the 
new predictions with those of the original Bau05 and Bow06 mod- 
els. To do this, we leave the parameter values for processes other 
than the SF law unchanged, to isolate the effects of changing the SF 
prescription. We apply the new SF laws consistently throughout the 
galaxy formation calculation (see Appendix[C]for an illustration of 
applying the SF laws to the original models over a single timestep). 

Fig. shows the evolution of the cosmic SFR density (top), 
distinguishing between quiescent (middle) and burst (bottom) SF 
modes, for the Bau05 (left) and the Bow06 (right) models. In the 
Bau05 model, the peak in SF activity is at z as 2.5. The peak lo- 
cation (though not its height) is somewhat sensitive to the choice 
of SF law, moving to higher redshift in the variant models. The ex- 
treme cases are the KS and KS. thresh SF laws, in which the peak 
shifts to z ~ 3.5. In the case of the KMT and BR.nonlin laws, the 
peak moves slightly to z ~ 2.8, while the BR law gives a very sim- 
ilar SFR density evolution to the original Bau05 model. The change 
in the location of the SFR density peak is driven by the change in 
the level of quiescent SF activity. SF laws that produce higher SFR 
densities in the quiescent mode at high redshift (i.e. KS, KS. thresh, 
KMT, BR.nonlin) are characterised by a SFR-gas density correla- 
tion Esfr — S gas steeper than linear (e.g. TV > 1 in Eq. |4j- A 
consequence of the larger quiescent SFR at high redshifts is less 



cold gas available to fuel starbursts, and so the SFR density in the 
burst mode drops. 

In the Bow06 model (right panels in Fig. Q3, the evolution of 
the total SFR density is largely unaffected by the choice of SF law 
even though the relative importance of the quiescent and burst SF 
modes changes by more than one order of magnitude. This is due to 
the inclusion of starbursts triggered by disk instabilities (§2). The 
smaller SFR density in the quiescent mode relative to the original 
model when using the BR, BR.nonlin and KMT SF laws (middle 
panel), leads to significant amounts of cold gas in galaxy disks. This 
results in a higher frequency of unstable disks and gas rich merg- 
ers, and also more prominent SBs due to the larger gas reservoirs. 
This increase in starburst activity compensates for the lower quies- 
cent activity. In the case of the KS and KS .thresh laws, the slightly 
smaller SFR density (by « — 0.1 dex) in the quiescent mode is not 
offset by an increase in the SF activity from disk instabilities. This 
is because the condition for a disk to become dynamically unstable 

— 1/2 

is only weakly dependent on mass (oc M diB ^ ). Hence a substan- 
tial change in disk mass is needed to make a galaxy unstable, and 
on average the change in disk mass is less than 0.1 dex for the KS 
and the KS. thresh laws, while for the other SF laws it is at least 
0.4 dex. Remarkably, all of the new SF laws produce peaks of burst 
and quiescent SF activity separated by at least Az w 1, while in 
the original Bow06 model, the peaks overlap. 

The differences on applying the new SF laws in the Bau05 and 
Bow06 models indicate the role of physical processes other than 
quiescent SF in shaping the evolution of galaxies, e.g. disk instabil- 
ities, and support our decision to leave the other model parameters 
unchanged. However, there are remarkable similarities between the 
predictions which result on using different SF laws. Even though 
the total density of SF activity is unchanged, the relative importance 
of the burst and quiescent modes depends sensitively on the choice 
of SF law, due mainly to the change in the cold gas content of galax- 
ies (see §4). The insensitivity of the total SFR d ensity to the choice 
of SF law has been noted in previous work (e. g. lStringer & Benson! 
120071; button et alj|20ld : ISchave et alj|20ich . This can be under- 
stood as resulting from the SFR (combined with the rate of gas 
ejection from galaxies due to SN feedback, which is proportional 
to the SFR) on average adjusting to balance the rate of accretion 
of cooling gas. The SFR will tend to adjust in this way whenever 
the gas consumption timescale tsf = M co \a/ SFR is less than the 
age of the universe at that redshift. A direct consequence of this 
interplay between the modes of SF is the insensitivity of the optical 
and near-IR galaxy LF and optical colours to the choice of SF law 
(see Appendix IdI for a comparison of the predicted LFs and other 
galaxy properties when using different SF laws). 



4 COLD GAS MASS CONTENT OF GALAXIES 

The cold gas mass content of galaxies is sensitive to the choice of 
SF law, since this determines the rate at which gas is converted into 
stars. Here we focus on three separate observational probes of gas: 
the cold gas mass function (CMF), the evolution of the cold gas 
density of the universe and the gas-to-luminosity ratio in galaxies. 

4.1 Cold gas mass function 

Fig. [2] shows the CMF at different redshifts for the original mod- 
els and with differ ent SF laws. The observed 2 = CMF of 
IZwaan et al.l d2005l) is plotted to highlight evolution. Note that the 
observed CMF plotted here is derived from the observed HI mass 
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Figure 1. The evolution of the cosmic star formation rate per unit volume. The left-hand panels show the Bau()5 model and its variants, and the right-hand 
panels show the Bow06 model. The top row shows the total star formation rate, the middle panels show the quiescent star formation and the bottom panels the 
star formation in bursts. For reference, in the middle and bottom panels the total SFR density in the original models is shown by a thick gray solid line. 



function assuming a constant H2/HI ratio of 0.4 and a hydro- 
gen mass fraction X — 0.7 4 (§2.3.1) to convert the HI measure- 
ments into cold gas masses ( Bau gh et alj|2004l : iPower et alj|201ol ; 
iKim et al.ll2010h . 

The Bau05 model is in quite good agreement with the ob- 
served CMF, except at the lowest gas masses. The parameters in the 
original model were chosen to match the gas-to-luminosity ratios 
discussed in ^4.31 The modified SF laws all lead to poorer agree- 
ment with observations. On the other hand, the original Bow06 
model is in poor agreement with the observed CMF. Some of the 
variant SF laws in this case lead to better agreement with the ob- 
served CMF. 

Most of the changes in the CMF with the new SF laws can be 
understood as resulting from the changed SF timescales in disks, 
defined as tsf = M co \d/ SFR. (Note that tsf would be equal to 



the gas depletion timescale if there was no accretion of fresh gas 
from cooling, no ejection of gas by SN feedback and no recycling 
of gas to the ISM by dying stars.) Since the treatment of gas cooling 
in halos is unchanged, accretion rates of cold gas onto disks are also 
essentially unchanged, so an increase in the SF timescale tsf leads 
to larger cold gas masses in disks and vice versa. The new SF laws 
depend on galaxy properties such as the radius, stellar mass and gas 
metallicity in different ways from the old laws. Furthermore, all of 
the new SF laws are non-linear in cold gas mass, at least in the low 
gas density regime, so that rg f always increases if M co td becomes 
low enough, in contrast to the original SF laws, for which tsf is 
independent of the remaining gas mass. As a result, the changes in 
gas content due to changing the SF law themselves depend on the 
galaxy mass and redshift (see Appendix [Cl for the change in SFR 
with galaxy mass and redshift for different SF laws). 
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Figure 2. The cold gas mass functions at z = and z = 2 for the Bau()5 
(upper two panels) and the Bow06 (lower two panels) models and for the 
different forms of the SF law. The line colours and styles are as in Fig.[T] For 
referen ce, the observational estimate of the z = CMF from lZwaan et alj 
12005I) is shown using grey symbols in all panels. A constant H2 /HI ratio 
of 0.4 has been assumed to convert the HI observations into cold gas masses. 



Starbursts triggered by disk instabilities have a bearing on gas 
content in the Bow06 model. Changing the SF law can lead to 
changes in the total disk mass, making the disk more prone to in- 
stabilities if the total mass is larger. Following an instability, all of 
the gas is consumed in a starburst. Thus, paradoxically, longer SF 
timescales in disks can lead to lower final gas contents in some 
cases. Starbursts in the Bau05 model are triggered only by galaxy 
mergers, and the frequency of these is essentially unaffected by 
changes in the gas content of galaxies. 

In the Bau05 model, the SF timescales are in most cases 
shorter with the new SF laws, and so the CMF is also lower at 
most gas masses. For the Bow06 model, the results of changing the 
SF law are more mixed, although the CMF is generally lower at 
intermediate gas masses. An interesting effect appears at low gas 
masses, where in all cases, the number density is larger than with 
the original SF laws, bringing the models closer to the observa- 
tional data at z — 0. This results from the new SF laws being non- 
linear in the cold gas mass, so that the gas is depleted less rapidly 
when the mass becomes very small than with the original (linear) 
SF laws. The effects of halo mass resolution are more severe for 
the Bow06 model, which uses N-body merger trees from the Mil- 
lennium simulation. This is what causes the turndown in the CMF 
forM cold < 1O 9 /i _2 M . 

The high-mass end of the CMF in the Bow06 model evolves 
significantly from z = 2 to z — due to the longer SF timescales 
at lower redshifts (see Eq. 2). With the new SF laws, this evolution 
is weaker, due in part to the higher frequency of disk instabilities 
triggering starbursts which consume the gas. 

The evolution of the CMF is therefore a useful diagnostic to 
distinguish between the SF laws, since they lead to changes that are 
differential with galaxy mass. 

4.2 Global cold gas density evolution 

The evolution of the global cold gas density, p co id, is shaped by 
the same processes as the CMF. Fig.[3]shows Pcoid as a function of 
redshift for the Bau05 (top panel) and Bow06 (bottom panel) mod- 
els and the variants with new SF laws. For reference, grey symbols 
show different observations as listed in the legends. In general, all 
the new SF laws predict lower p co id than the original models across 
the whole redshift range. 

In the Bau05 model, the reduced p co id with the new SF laws is 
due to the shorter quiescent SF timescales in low and intermediate 
cold gas mass objects (see right-hand panels of Fig. ICU . In gen- 
eral, the shorter the quiescent SF timescale (on average) the smaller 
Pcoid compared to the original model. Note the jump of p co id at 
z ~ 6 in the Bau05 model is due to the assumed reionization red- 
shift Zreion = 6 (§6.1). 

On the other hand, in the Bow06 model the BR, BR.nonlin 
and the KMT SF laws at very high redshift (i.e. 2^6) give similar 
p co id to the original model, mainly due to the compensation be- 
tween quiescent and burst SF activity. The BR, BR.nonlin and the 
KMT SF laws give reasonable agreement with the observed p co id 
at z = 0, improving over the original model. 

The Bow06 model shows a more rapid decrease of p co id with 
increasing redshift than the Bau05 model. This results in part from 
the stronger SNe feedback in the Bow06 model, but is also due to 
the difference in the mass resolution of the halo merger trees used 
in the two models (see ^2. H . The dotted line in the bottom panel 
of Fig. [3] shows the effect on the standard Bow06 model of using 
Monte Carlo trees with a minimum mass lO lo /i- 1 M (l + z)" 3 
instead of N-body trees with a fixed mass resolution (shown by the 
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Figure 3. The cold gas mass density in the universe, p co id . in units of the 
critical density at z = 0, p C rit,z=o> as a function of redshift. The results 
for the original Bau05 and Bow06 models and with the new SF laws are 
shown in the top and bo ttom panels r espectively. Grey symbols show ob- 
servational results from Zwaan et alj [ 19971 120051 op en and filled circles, 
respectively), [Rao & Briggsl 1 19931 open triangles) and lPeroux et alj fc003l 
open and filled squares). For the latter, open symbols indicate the cold gas 
density inferred from damped Lyo systems, while filled symbols include 
a correction for gas clouds with lower column density than originally de- 
tected. The dotted line in the bottom panel shows the effect of changing the 
halo mass resolution from the default values, as discussed in the text. 



Figure 4. Hydrogen gas mass to luminosity ratios, Mh/Lb, as a function 
of ZJ-band magnitude for late-type galaxies. Late-type galaxies in the model 
are those with a bulge-to-total luminosity in the B-band of B/T ^ 0.4. 
The Bau05 and Bow06 models are shown in the top and bottom panels 
respectively. The line colours and styles which represent the different SF 
laws are the same as in Fig. ID II For cl arity we only show the media n of 
each model . Observational results from Huchtmeier & Richtei] Jl988l) and 
ISaed (7993) are shown as small gray symbols. The triangles with errorbars 
show the median and the f and 90 percentiles for the observations. In order 
to compare directly with the observations, we assume that 74% of the cold 
gas mass predicted by the models is hydrogen. 



solid black line). The effect of improving the halo mass resolution 
is larger at high redshifts. 



4.3 Gas-to-luminosity ratios of galaxies 

Fig. [4] shows the gas-to-luminosity ratio as a function of absolute 
magnitude in the B-band at z = for late-type galaxies (those with 
bulge-to-total luminosity in the B-band, B/T ^ 0.4 in order to 
make a fair comparison with observations), for the original models 
and with the new SF laws. 

The Bau05 model predicts gas-to-luminosity ratios in good 
agreement with the observations, as these data were used as the 
main constraint on the parameters in the original SF recipe (eqn(3]l- 
The new SF laws in the Bau05 model predict gas-to-luminosity ra- 
tios which increase with luminosity, contrary to the observed ten- 



dency. This is due to the new SF laws predicting generally shorter 
quiescent SF timescales in low and intermediate mass galaxies but 
longer timescales in high mass galaxies, when compared to the 
original model (see Appendix |Cl>, 

The original Bow06 model predicts gas-to-luminosity ratios 
which are in poor agreement with the observations, except at the 
highest luminosities. However, the new SF laws generally pre- 
dict lower gas-to-luminosity ratios for galaxies fainter than Mb — 
51og(/i) » —20, due to shorter quiescent SF timescales, bringing 
them into better agreement with the observations. 
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5 THE EVOLUTION OF GALAXIES IN THE SFR VS. 
STELLAR MASS PLANE 

The level of SF activity in galaxies at any redshift is dependent on 
the available cold gas and the timescale on which the gas is con- 
sumed (which is set by the choice of SF law). Here we examine the 
influence of different SF laws on the distribution of galaxies in the 
SFR-stellar mass (SFR- A/*) plane and its evolution with redshift. 
The distribution of galaxies in this plane offers a way to constrain 

the form of the SF law. 

Brinchman n et al.l fc004l) found that local star-forming galax- 
ies lie on a relati v ely narrow r egion in t he SFR-AL plane . 
iNoeske et al] fe007t) . lElbaz et alj J2007I) and iDaddi et alj d2007h 
extended this work to higher redshifts, finding that the "main se- 
quence" of star-forming galaxies in the SFR- A/* plane remains in 
place up to at least z ~ 2 with a similar slo pe to that at z ~ . 
The zero-point evolves strongly with redshift. [Santini et al.1 ( 2009) 
found evidence for a bimodal distribution of galaxies in the SFR- 
M* plane, with a "passive" sequence below that of actively star- 
forming galaxies. Below, we investigate whether or not the models 
are able to explain the observations, and to what extent the distri- 
bution of galaxies in the SFR- A/* constrains the SF law. 



5.1 The local star formation rate-stellar mass plane 

Fig. [5] shows the predicted distribution of galaxies in the SFR-M* 
plane at z = for the Bau05 and Bow06 models and for variants 
using the new SF laws. The colour shading indicates the regions 
within which 40% to 99% of the model galaxies fall in the SFR-M* 
plane, when the distributions are volume-weighted and normalized 
in bins of stellar mass. We also show, using solid and dotted black 
contours, the regions of the plane where satellite and burst galax- 
ies respectively constitute more than 50% of the population within 
cells in the SFR- A/* plane. Here we define burst galaxies as those 
in which the SB mode accounts for at least 20% of their total SFR. 
Most of the galaxies (> 99% when integrated over the whole SFR- 
Mi, plane) are undergoing quiescent SF. Some of these contours 
break up into "islands", but this is an artifact of some cells in the 
SFR-A/* plane containing a small number of galaxies. 

We see from Fig.|5]that all of the models apart from the Bow06 
model show two sequences in the SFR-M* plane, although the 
shapes and widths of these sequences vary significantly depend- 
ing on the SF law. We will call the upper sequence with higher 
SFRs at a given M* the "active" or "star-forming" sequence, and 
the lower sequence with lower SFRs the "passive" sequence. In the 
original Bow06 model there appears to be only an active sequence, 
but there is a large spread of galaxies down to lower SFRs. The 
two sequences appear more distinct for the variants of the Bow06 
model with the new SF laws than for the Bau05 model and its vari- 
ants. For the Bow06 model and variants, the active SF sequence ex- 
tends down to the lowest stellar masses plotted, while for the Bau05 
model and its variants, the active sequence is only present for higher 
stellar masses. We now examine these differences in more detail. 

In the Bau05 model and its variants, there are comparable 
numbers of active and passive galaxies for AL > 10 10 /i _1 M o . 
At lower masses most of the galaxies are on the passive sequence, 
with the active sequence completely disappearing for M+ < 
lO 9 /i -1 Af0. The passive sequence is generally broader than the ac- 
tive one. The slope of the passive sequence is generally steeper with 
the new SF laws, producing lower SFRs at M* < 10 10 h' 1 M Q . 
This is closely related to the lower gas content in these galaxies, 
as seen in Fig. [4] The exception is the KS .thresh law, where the 
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Figure 6. SFR vs. stellar mass plane for galaxies at z = in the Bow06 
model with the different SF laws tested. The contours outline the 95% per- 
centile regions for satellite (red) and central (blue) galaxies, after normaliz- 
ing in bins of stellar mass. 



passive sequence disappears at low SFR and low M*. This results 
from the cut-off in SF for galaxies with gas surface densities be- 
low the critical density (see Eq. |5j- Indeed, 70% of all galaxies in 
this model with M* > 10 7 /i _1 M have SFR = 0. However, in 
galaxies with Af* > 10 9 ft _1 M , this percentage falls to 30%. 
This means that the critical density cut-off completely switches off 
SF activity in most low mass galaxies at z = 0. 

For the variants of the Bow06 model the passive sequence 
appears generally narrower and more offset from the active one, 
compared to the Bau05 model variants; furthermore, there seems 
to be more of dependence of the shape on the SF law. As for 
the Bau05.KS. thresh model, the passive sequence disappears com- 
pletely at low SFRs in the Bow06.KS. thresh model. 

In order to gain further insight into the nature of the upper and 
lower sequences in the SFR- A/* plane, we split the sample into cen- 
tral and satellite galaxies. Fig.[6]shows the 95% percentile contours 
after normalization in bins of stellar mass, for central (blue con- 
tours) and satellite galaxies (red contours) for the Bow06 model 
and its variants. Satellite galaxies are found predominantly on the 
lower sequence, while the upper sequence is dominated by central 
galaxies, as expected from the solid black contours in Fig. [5] This 
implies that the inflow of newly cooled gas onto central galaxies 
is a key process in shaping the "active"sequence, while the rate 
at which the cold gas reservoirs in satellite galaxies are consumed 
shapes the lower, "passive" sequence. In satellite galaxies, the new 
SF laws applied to the Bow06 model typically produce longer SF 
timescales (see Appendix|C]l. This is a consequence of the new SF 
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Figure 5. Distribution of galaxies in the SFR vs. stellar mass plane at z = for the Bau()5 and the Bow06 models and their variants with the new SF laws 
(left and right-hand sets of panels respectively). Each panel corresponds to a different model, as labelled. The coloured contours show the regions within which 
different fractions of the galaxies lie for given stellar mass, volume-weighted and normalized in bins of stellar mass, with the scale shown by the key. We show, 
using dotted and solid black contours, the regions in which starburst and satellite galaxies respectively make up more than 50% of the galaxy population. 



laws being non-linear in cold gas mass, hence the SF timescale 
becomes long when the gas mass becomes small, which happens 
when galaxies no longer accrete newly cooled gas. Therefore satel- 
lite galaxies consume their gas reservoirs more slowly than in the 
original model, which results in larger cold gas masses and higher 
SFRs at later epochs, and allows the galaxies to remain on a well 
defined sequence for longer. In contrast, with the old SF recipe the 
gas reservoirs of satellites decline exponentially with time. 

The change in the slope of the passive sequence at 
~ W 10 h~ 1 M El in the Bow06.BR, Bow06.BR.nonlin and 
Bow06.KMT models is due to the change in the dominant regime 
in the Esfr, — E gas relation. Low mass galaxies are predominantly 
in the low E gas regime, for which most of the hydrogen is in the 
atomic phase, and the SFR has a steep dependence on E gas . Indeed, 
the steeper the relation in the Esfr — E gas plane at low E gas , the 
flatter the passive sequence at low stellar masses (see Fig.|5]for the 
KMT SF law). However, this simple picture only holds for satellite 
galaxies, since centrals have a constant supply of newly cooled gas 
which shapes the upper sequence in the SFR- A/* plane, making it 
insensitive to the exact choice of SF law. 

In the Bau05 model the satellite and central galaxy sequences 
overlap more due to the longer timescale for gas to be reincorpo- 
rated into the host halo after ejection by SNe, compared with the 
Bow06 model (see §2). This is particularly important at low stel- 
lar masses, for which SNe feedback is very efficient at ejecting gas 
which has previously cooled. The long reincorporation timescale 
means that star formation is nearly shut down in most low-mass 



central galaxies. Hence the active SF sequence disappears at low 
M^. The passive sequence at low M* therefore includes both satel- 
lite and central galaxies. We confirmed this by carrying out the ex- 
ercise of applying the Bow06 prescription for the gas reincorpora- 
tion timescale in the Bau05 model, and found an active sequence 
which looks similar to that in the Bow06 model. The slope and dis- 
persion of the active SF sequence are insensitive to the details of 
the SNe feedback (or the SF law; see Fig. [5), since making feed- 
back stronger or weaker does not change its shape. We conclude 
that the reincorporation timescale of the ejected gas is the main 
process controlling the form of the active SF sequence, while the 
passive sequence is mostly dependent on the choice of SF law. This 
relation offers a new way to constrain the SF law. 



5.2 Comparison with observations of the SFR-M, plane at 

2 = 

Indications of multiple sequence s in the SFR-M* p l ane h ave been 
found in the local Universe by iBrinchmann et al. fc004l). a nd at 
high redshifts by IPamen et al.l d2009h . ISantini et al.l j2009h and 
Rodig hiero" et alj < l2010h " However, all of these studies are affected 
to some extent by selection effects against the inclusion of galax- 
ies with low SFRs and/or low stellar masses, and these must be 
accounted for in any comparison with models. We start by compar- 
ing with observations of local galaxies in this subsection, and then 
consider higher redshifts in the next subsection. 

In the top panel of Fig.[7]we show the observed distribution of 



The impact of star formation laws 1 1 



99% 95% 90% 80% 64% 50% 




8 9 10 11 

iog(M stellar /h-' ivy 



12 



2 

1 



-1 

-2 
-3 

2 
1 

~ 

-1 

©-2 

I" 3 
2 

"35 1 
* 
-1 
-2 
-3 

2 
1 

-1 
-2 
-3 



\ Bau05 J 

f^M''V!| | | | | 


f- Bow06 i 

;, ,,,-mVm i i i ; 


\ Bau05.KMT ^^g^^ 


f Bow06.KMT 1 


j. .".Vi 1 1 1 j 

r Bau05.KS.thresh \ 


^..'iVi| 1 1 1 j 

r Bow06.KS.thresh i 

; -mVm| i i i ; 


; i i i i 

r^Bau05^R^^^^^^^^ 


r Bow06.BR -) 



8 9 

iog(M stellar 



10 11 

/ h ' Mq) 



8 9 

iog(M stellar 



10 11 

/ h 1 ivy 



Figure 7. Top panel: Observed distribution of galaxies in the SFR vs. stel- 
lar mass plane for the SPSS DR7 spectrosc opic sample of galaxies, updat- 
ing the analysis by Brinchmann et al. 1 2004). The distributions are volume- 
weighted and normalized in bins of stellar mass as in Fig. [5] and the shading 
shows the regions within which different fractions of galaxies lie at each 
stellar mass. The black line shows the effective SFR sensitivity, which de- 
pends on stellar mass (Brinchmann, priv. com.). Bottom panel The same as 
the top panel for model g alaxies at z = 0.1 after applying the cuts used by 
iBrinchmann et alj (2004) in the SDSS DR7 sample, for different models as 
labelled in each panel. For reference, the dotted lines show the sensitivity 
limit of the SFR estimates in the SDSS DR7. The shading is as in the top 
panel. The KS and BR.nonlin models are not shown due to their similarity 
with the KS. thresh and BR models respectively. 



galaxies in the SFR-Af* plane based on the SD SS Data Release 7 
(DR7f[ which corresponds to an update of the Brinchmann et al. 
1 20041) analysis. Stellar masses are determined from spectra and 



4 Data was downloaded from the public 

http : / /www . mpa-garching . mpg . de/ SDSS/DR7/. 



webpage 



broad band magnitudes following Kauffma nn et alj d2003l) . while 
SFRs are derived primar ily from the H a emission line following 
IBrinchmann et al.l d2004h . We refer the reader to these papers for 
further details. We have constructed contours of the distribution 
for the subsample of galaxies in the spectroscopic catalogue which 
have estimated SFRs, after volume-weighting and normalizing in 
bins of stellar mass. As in Fig. [5] the colour shading indicates the 
regions within which different fractions of galaxies lie for a given 
stellar mass. The solid black line shows the approximate SFR sen- 
sitivity limit in the SDSS analysis (Brinchmann, priv. comm.). The 
minimum measurable SFR depends on stellar mass both because 
it depends on detecting spectral features above the stellar contin- 
uum, and because galaxies of higher M* in the SDSS sample tend 
t o lie at la r ger di stances. Since the SDSS estimat es are based on 
a lChabrierl (2003) IMF for the stellar masses, and a lKroupal ( l200ll) 
IMF for the SFRs, we rescale the stellar masses in t he SDSS sample 
by a factor of 0.89 and the SF Rs by a fac t or of 1 .1 dBell et alj20"03h 
to correspond to the use of a lKennicutu dl983l) IMF for quiescent 
star formation in the models. (In the Bau05 model starbursts form 
stars with a top-heavy IMF, but this only affects a small fraction of 
the galaxies in the SFR-M* plots. SB galaxies are typically located 
above the colour contours at higher SFRs, see dotted lines in Fig. [5] 
and Fig. [8}. 

To make the comparison with the SDSS data simpler and 
fairer, we plot the model predictions in the SFR-Ai* plane again 
in the bottom panel of Fig [7j but this time imposing some addi- 
tional selections to better match the SDSS sample. Since the me- 
dian redshift of the SDSS spectroscopic sample is around z = 0.1, 
we create a sample of model galaxies for z = 0.1, and then apply 
the apparent magnitude limit r < 17.77 of the SDSS spectroscopic 
sample. We also apply the SFR sensitivity limit of the SDSS sample 
shown in the top panel of Fig. [7] Fig UJ shows the SFR- M* plane 
for the Bau05 (left panels) and the Bow06 (right panels) models 
at 2 = 0.1 after applying these cuts. Note that we do not show 
the BR.nonlin and KS models given their similarity to the BR and 
KS. thresh, respectively. 

For the SDSS sample plotted in Fig. [7] we see that most galax- 
ies in the mass range 10 7 < Af* < 10 10 h~ x Mq lie on a single 
star-forming sequence. At high masses > 10 10 /i _1 A/©, a sec- 
ond sequence appears, roughly parallel to the first but offset below 
it. Comparing with the model distributions in Fig. UJ the general 
distribution of galaxies in the SFR-Af* plane seems better matched 
by the variants of the Bow06 model rather than by any of the vari- 
ants of the Bau05 model. All of the variants of the Bow06 model 
produce an upper star-forming sequence covering the whole mass 
range 10 7 < M* < 10 11 h~ 1 M Q , with similar slope to the ob- 
served one, though with a smaller dispersion, particularly at low 
masses. The original Bow06 model does not reproduce the second 
sequence seen in the SDSS at high masses, but most of the vari- 
ants with new SF laws do produce a similar feature. On the other 
hand, in the Bau05 model, the distribution of galaxies in the SFR- 
Mi, plane appears rather different from the observed one. All the 
variants show an upper sequence of star-forming galaxies, which 
has a similar slope and amplitude to that seen in SDSS. However, 
this does not extend to low masses, unlike in the SDSS. The most 
prominent feature for the Bau05 model and its variants is the broad 
"passive" sequence lying below the "active" star-forming sequence, 
and extending down to the lowest masses plotted. This does not cor- 
respond to what is seen in the SDSS sample, where all low-mass 
galaxies appear to lie close to the "active" star-forming sequence. 
The main origin of these differences between the Bau05 and Bow06 
model for low stellar masses is the difference in the timescale for 
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Figure 8. Distribution of galaxies in the SFR vs stellar mass plane at z = 1 
(left column), z = 3 (middle column) and z = 6 (right column) for the 
Bau05 (top panel), the Bau05.BR (second panel), the Bow()6 (third panel) 
and the Bow06.BR models (bottom panel). The shading shows the distribu- 
tion of galaxies in this plane, normalized in bins of stellar mass. The thick 
da shed straight lines show the obser vations of the "SF s equences" reported 
bv lElbazetal] 120071) at z ~ 1 and IStark et"aj] l2009j) at z ~ 6, and are 
plotted over the mass ranges probed by the observations. We also show re- 
gions in which starburst and satellite galaxies make up more than 50% of 
the population using dotted and solid black contours, respectively. 



the reincorporation into halos of gas which has been ejected by su- 
pernovae (see 6.1). 

The passive sequence seen in Fig. [5] for the variants of the 
Bow06 model falls below the SFR sensitivity limit of the SDSS 
study at lower masses, A/* < lO lo ft -1 M0. However, observa- 
tional indications of a passive gala xy sequence at low masses have 
been reported bv lWoo et al.l d2008h for local group dwarf galaxies. 
These objects follow a distinct SF sequence with much lower SFRs 
than the main star-forming sequence seen in the SDSS sample, sim- 
ilar to our predictions for passive gala xies in the Bow06 mo del 
variants with the new SF laws. Similarlv. lSkillman et alj ( 120030 re- 
ported SFRs and stellar masses for galaxies in the Coma cluster 
which agree with the lower sequence seen in the top panel of Fig.[7] 
and also overlap with our prediction for massive satellite galax- 
ies in the Bow06.BR, Bow06.BR.nonlin and Bow06.KMT models 
(Fig.©. 



5.3 The SFR-M* relation at high redshift 

Evidence for a star-forming sequence in the SFR- Af* plane has also 
been found at high redshifts. The slope of the observed relation 



seems to be constant over the redshift range prob ed, while the zero 
point evolves strongly from z ~ to z ~ 4 ( lElbaz et alJliooH : 
iDaddi et alj2007l :|p erez-Gonzalez et al. 2008), but appears roughly 
constant for 4 < z < 6 JStark et alj |2009). However, we caution 
that the selection effects in some of the high-redshift samples are 
very strong, since some are selected on the basis of SFRs. Further- 
more, the observational estimates of SFR and stellar mass typically 
become more uncertain at high er-z, due to the mor e limited data 
available and the effects of dust JStringeretaI1l201lh . 

Both models in all of their variants predict an active SF se- 
quence in good agreement with observations. The features of the 
models in the SFR- M* plane at z — (see Fig. [5]l are preserved 
up to z « 2. At z > 2, the passive sequence starts to disappear, 
implying that at higher redshifts most of the galaxy population is 
undergoing vigorous SF and lies on the active SF sequence. To il- 
lustrate this general behaviour, Fig. [8] shows the predicted distri- 
bution of galaxies in the SFR- A/* at three different redshifts for 
the Bau05 and the Bow06 models, together with the BR variant. 
Also shown are the medians of the observed relations at different 
redshi fts, plotted as thick dashed linefl 

ISantini et all J2009h found that bimodality in the SFR-M* 
plane is clearly present in massive galaxies at z < 1, but becomes 
weaker as redshift increases, being almost absent at z ~ 2. This 
appears consistent with the predictions of the Bow06 model with 
the new SF laws, but is probably less so with predictions from the 
Bau05 model. However, the SFR sensitivity limit in the sample of 
Santini et al. is quite high (SFR > 1O~ 2 M yr -1 ) relative to the 
model predictions at lower masses, making it difficult to study the 
passive population at high redshift. 

The new SF laws, through their impact on the cold gas con- 
tents of galaxies, strongly affect the passive population and its evo- 
lution as seen in the SFR-A/ + plane, in a way that could be used 
as an important constraint on the SF model, if observational data 
could probe low enough SFRs. On the other hand, the active SF 
sequence is much less affected by the choice of SF law, but instead 
depends on the assumption about how quickly gas ejected by SN 
feedback is reincorporated into the host halo. 



6 DISCUSSION AND CONCLUSIONS 

Improvements in the quality and quantity of observations of the 
spatially-resolved star formation rates an d atomic and molec- 
ular gas content s in local galaxies (e.g. iKennicutt et alj 120071 : 
Bigiel et al. 2008) have allowed the development o f both empirical 
(e.g. lBlitz & Rosolowsk\]|2006l:lLerov et alJ|2008h and theoretical 
star formation (SF) models (Kru mholz et all2009bf) . These models 
prescribe a relation between the surface density of star formation 
Esfr and gas S gas , often depending also on other physical quanti- 
ties such as the stellar surface density and gas metallicity. We have 
revisited the SF recipes for quiescent galaxy formation in galaxy 
disks used in the galaxy formation model GALFORM and imple- 
ment new, parameter-free SF laws with the aim of isolating obser- 
vations that could be used to distinguish between SF models (see 
lLagos et al. 201 1 for the predictions of these models for the atomic 
and molecular hydrogen content of galaxies). 

We test the following SF laws: (i) the Schmidt law, which has 

5 Since in the observational samples a universal Salpeter IMF was as- 
sumed when estimating SFRs and stellar masses, we scale thei r stellar 
masse s and SFRs by fac tors of 0.5 and .91 respectively (e.g. iBell et alj 
120031) . to correspond to a lKennicutJ Jl983l) IMF. 
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the form E SF r oc E^ s with N = 1.4 (the KS law) and (ii) its vari- 
ant, the Kennicutt-Schmidt law, which includes a SF cut-off at low 
gas su rface densities , motivated by gravitational stability consider- 
ations jToomrel 19641 ; the KS. thresh law); (iii) A SF law of the form 
Esfr oc ^sf /moi S gaa , where / mo i depends on the hydrostatic ga s 
pressure (see Eq. Blitz & Rosolowskvl l2006l ; lLerov et al.1 12008). 
and z^sf is an inverse timescale of SF activity in molecular gas, 
assumed to be constant (the BR law); (iv) A variant of the previ- 
ous SF law, but assuming ^sf = « / SF(S gas ) which depends on the 
total gas surface dens ity (see Eq. |9l the BR .nonlin law), (v) The 
theoretical SF law of lKrumholz et al.l (2009b) in which the ratio of 
surface densities of molecular to atomic gas depends on the bal- 
ance between the dissociation of molecules due to the interstellar 
far-UV radiation, and their formation on the surface of dust grains 
(the KMT law). We apply the new SF laws to the quiescent SF 
mode which takes place in galactic disks (see §2.1). 

We have applied the new SF laws to two variants of 
GALFORM, those of Baugh et al. (2005; Bau05) and Bower et al. 
(2006; Bow06). These models have many differences in their in- 
put physics, including different laws for quiescent SF, and further- 
more their parameters were tuned to reproduce different sets of ob- 
servational constraints. We have left the parameters for all other 
physical processes, including SF in starbursts and feedback from 
supernovae and AGN, unchanged when running with the new qui- 
escent SF laws, rather than retuning these other parameters to try to 
reproduce the original observational constraints. This allows us to 
isolate the impact of invoking different SF laws. Previous attempts 
to include similar SF recipes in semi-analytical models have not 
focused on how galaxy properties are affected by changing the SF 
law al one, but have also made changes to other model parameters 
(e.g. lDutton et al.ll201Cl ; lc~ook et alj20ld : lFu et alfeOld) . 

The choice of quiescent SF law by construction affects the 
quiescent SF activity. In the most extreme case, the cosmic SFR 
density in the quiescent mode changes by an order of magnitude 
when compared with the original model. However, the total SFR 
density evolution is remarkably insensitive to the choice of SF law. 
In the case of the Bau05 model the new SF laws increase the SFR 
in the quiescent mode, while the contrary happens in the Bow06 
model. This is compensated for by weaker starburst activity in the 
Bau05 model and more vigorous starbursts in the Bow06 model. 
This adjustment is particularly effective in the case of the Bow06 
model, which includes starbursts triggered by disk instabilities. A 
related consequence is the insensitivity of the present-day bj-band 
LF and g — r colour distributions and the evolution of the A'-band 
LF to the choice of quiescent SF law. 

Since the gas is consumed at different rates during quiescent 
and starburst SF, the change in the quiescent SF law leads to large 
differences in the cold gas content of galaxies. For both models the 
new SF laws predict more galaxies with low cold gas masses, fewer 
with intermediate cold gas masses and similar or fewer with high 
cold masses. This leads to a lower cosmic mean cold gas density at 
all redshifts. For the Bow06 model, these chang es help to improv e 
the agreement with the local observations bv lZwaan et alj J2005h . 
We find similar changes in the Mh/Lb gas-to-light ratios as func- 
tions of luminosity. 

Finally, we investigated the distribution of galaxies in the 
SFR-A/* plane at different redshifts for both the old and new SF 
laws. In most cases, the models predict two sequences in this plane, 
an "active" sequence with higher SFRs and a "passive" sequence 
with lower SFRs. The form of the "active" sequence, which is dom- 
inated by central galaxies, appears insensitive to the quiescent SF 
law. This can be understood in terms of a rough balance being set 



up between accretion of gas due to cooling in the halo and con- 
sumption of gas by star formation and ejection by supernova feed- 
back (at a rate which is proportional to the SFR). The cold gas 
mass in a galaxy therefore adjusts to achieve this balance (for a 
given SF law), while the SFR itself is insensitive to the SF law. The 
form of the active sequence is, however, sensitive to the timescale 
for gas which has been ejected by supernova feedback to be rein- 
corporated into the halo, particularly in low-mass galaxies where 
the supernova feedback is very effective at ejecting gas. This rein- 
corporation timescale is relatively short (around a halo dynamical 
time) in the Bow06 model, but longer in the Bau05 model. As a 
consequence, the active sequence extends down to very low stellar 
masses in the Bow06 model and its variants with different SF laws, 
but it disappears at low masses in the Bau05 model and its variants. 

The form of the "passive" sequence, where most satellite 
galaxies are found, is much more sensitive to the SF law. Satellite 
galaxies are assumed to not accrete any cold gas, thus they are sim- 
ply using up existing reservoirs of cold gas by star formation and 
feedback, and their gas contents and SFRs are thus sensitive to the 
rate at which this has happened over their history. This depends on 
the SF law and weakly on other parameters, such as SN feedback. 
We propose that the direct measurement of the (low) SFR of pas- 
sive galaxies offers a constraint on the SF law and can be directly 
compared with the predictions made in this work. 

We made a qualitative comparison of the predicted distribution 
of galaxies in the SFR- M* plane with observational d ata. For z ~ 
0, we compared with the SDSS DR7 results of Brinchman n et al.l 
d2004l) which show two sequences, of active and passive galax- 
ies, at high stellar masses. The Bow06 model variants (particularly 
those with the BR, BR.nonlin or the KMT SF laws) appear in better 
agreement with these data than the Bau05 model, after accounting 
for the sensitivity limits. Observations at higher redshifts indicate 
the presence of an acti ve SF sequence that evolves in normalization 
but not in slope (e.g. [Brinchmann et alj|2004l ; iDaddi et al.ll2007l ; 
lElbaz et alj|2007l ; lstark et alj|2009l) . The model predictions appear 
in general agreement with these obs ervations for both the B au05 
and Bow06 models and their variants. iRodighieroetalJ l feOlOl) and 
ISantini et alj J2009h have found indications of bimodality at z ~ 1 
in massive galaxies, consistent with our predictions for the Bow06 
model in the three variants already mentioned. 

Future observations of SFRs, stellar masses and particularly 
the gas contents of galaxies at high redshifts will be crucial for test- 
ing whether current ideas about star formation laws applied at ear- 
lier stages in galaxy evolution. Telescopes such as ALMA, JWST 
and ELTs will make possible measurements of very low SFRs such 
as those predicted here for the passive galaxy population, while 
observations with ALMA, ASKAP, MeerKAT and eventually SKA 
will be critical for measuring the molecular and atomic gas contents 
of high-redshift galaxies. In this way we will be able to confirm or 
rule out the predictions of the models presented here. 
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APPENDIX A: NUMERICAL INTEGRATION OF THE 
STAR FORMATION EQUATIONS 

Here we outline the numerical solution of the equations that de- 
scribe the changes in the baryonic content of galaxies, namely the 
hot gas, cold gas, stellar mass and the metals in each component. 



Al The instantaneous SFR 



The radial profiles of galact ic disks are not resolved in GALFORM 
(except in special cases e.g.,Stringer & Benson 2 007h. We assum e 
disks are described by an exponential profile (see Cole et alj2000h . 



Sdisk(R) = E e 



-R/r 



(Al) 



where r c g is the scale length of the disk, Eo = Afdisk/27rr^ ff is the 
central surface density and Mdisk is the disk mass in the component 
traced (i.e. the cold gas mass, M co id for E gas , or the stellar mass of 
the disk, M ste n ari disk for E stc n ar ). Note that r c ff in an exponential 
disk is related to the half mass radius by t$o = 1.67r c ff . 

The instantaneous SFR in a given episode will be the sur- 
face integral of Esfr over the full disk, except for the case of 
the KS. thresh SF law, in which we integrate only over the unstable 
region, which corresponds to the solutions of E gas = E cr i t (see 
Eq.[5]l. In the case that the integral cannot be solved analytically (as 
in the BR, BR.nonlin and KMT SF laws), we use Romberg integra- 
tion to compute the global SFR by integrating Espr(R) over the 
whole disk. This method uses adaptive, equally-spaced divisions in 
R to achieve the required accuracy in each calculation. For simplic- 
ity, we do not distribute the newly formed stars only over the part 
of the disk where the SF activity occurs, but instead we assume 
that gas an d stars are always distrib uted with the same exponen- 
tial profile. S tringer & Benson] d2007l) found this to be a reasonable 
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approximation for the gas component in their radially resolved cal- 
culations of galactic disk evolution. However, this approximation 
could overestimate the radial extent of the stellar component of the 
disk in gas dominated galaxies (see IStringer & B enson 20Q7|; we 
revisit this point in Appendix C2). 



A2 The SF equations 

The SF activity in GALFORM is regulated by three channels: (i) 
accretion of gas which cools from hot gas haloes onto the disk, (ii) 
SF from the cold gas and, (iii) reheating and ejection of gas due 
to SNe feedback. These channels modify the mass and metallicity 
of each of the mass components (i.e. stellar mass, A'h, cold gas 
mass, Mcotd, hot halo gas mass, Mhot, and their respective masses 
in metals, M? , M ^ n]d and M^ nt ). T he system of equations relating 



these quantities is dCole et 



dCole et alfe OOO): 



M cold 
M hot 

MflH 



Mi 



(1 - R)^ 

Afcooi - (1 - R + P)il> 

-M coo l + Ptj) 

(1 - R)Z cM 4> 

Afcool^hot 

(p - (1 + p - R)Z cold )iji 
-Mcooi Z hot + (pe + /3Z cold )ip. 



(A2) 
(A3) 
(A4) 
(A5) 

(A6) 
(A7) 



Here p denotes the yield (the fraction of mass converted into stars 
that is returned to the ISM in the form of metals), e is the fraction of 
newly produced metals ejected directly from the stellar disk to the 
hot gas phase, R is the fraction of mass recycled to the ISM (in the 
form of stellar winds and SN explosions), Z co \ d — M^ oid /M cold 
is the metallicity of the cold gas, and /3 is the efficiency of stel- 
lar fe^dtec^TJescrirjtions t he feedback mechan isms can be found 
in iBenson et al.l d2003l) and lLacev et alj {2008). These equations 
are defined as a function of the instantaneous SFR, ip, that in 
the standard version of GALFORM has the form given by Eq. [T] 
with ip oc M co id. This form of the SFR has the advantage that 
Eqs. 15 — 20 can be solved analytically assuming that Af coo i is 
constant over a timestep. 

The new SF laws are non-linear in A-/ co id, so they require 
Eqs. (A2)-(A7) to be solved numerically. Once ip (§2.3.1) is cal- 
culated, we proceed to integrate Eqs. (A2)-(A7). In GALFORM, the 
equations tracking the evolution of the baryons are integrated over 
"halo timesteps" given by the time resolution at which the halo 
merger tree is stored, which in the case of both N-body and Monte 
Carlo trees is independent of the individual SF timescales in each 
galaxy. Hence, for some cases, the halo timestep could be large 
compared to the SF timescale. This makes it necessary to solve 
Eqs. (A2)-(A7) using adaptive stepsizes to achieve accurate so- 
lutions. We use the fourth-order Run ge-Kutta (RK) method with 
adaptive stepsizes dPress et al] 1 19921) . The quantities involved in 
the calculation of the instantaneous SFR, such as stellar mass, gas 
mass, metallicity of the gas, are updated in each sub-step. However, 
for simplicity, we assume that the disk scale length r c g and M coo i 
(see Eq. lAlt . remain constant during the integration over each halo 
step. After integrating Eqs. 15 — 20, we infer the galaxy luminosity 
by integrating the SF history over a halo timestep, interpolating the 
SFR between the values output at the substeps of the RK integra- 
tion. 



APPENDIX B: THE STAR FORMATION LAWS 

In this appendix we explain more extensively some physical details 
of the new SF laws included in the GALFORM model. 



Bl The critical surface density of the K98 star formation law 

In a thin isothermal disk, the critical surface density for gravita- 
tional instability of axisymmetric perturbations, E cr it, is given by 



Ecrit — 



(Bl) 



Qcrit VrG' 

where Qcrit is a dimensionless constant ~ 1 dToomrell 19641) , cr g is 
the velocity dispersion of the gas and k is the epicyclic frequency 
of the disk. For realistic gas/stellar disks we expect Qcrit > 1 due 
to the effects of non-axisymmetric instabilities and the gravity of 
the stars (see references in K89). It has been shown that <r g is ap- 
proximately constant within disks. This means that the radial de- 
pendence of Ecrit is determined by that of k, and can be estimated 
directly from a galaxy rotation curve, 



V2 



1 + 



VdR 



1/2 



For a flat rotation curve, 



Ecrit 



V2 



V 



Qcrit 7TG * R' 



(B2) 



(B3) 



K89 compared E cr it with the gas profile of the disk (E gas ), 
and found that at the radius of the outermost HII regions (which 



indicate recent SF activity), E gas /E c 



1.9 - 3.3 for Q c 



1 (after scaling by our choice of <r g ). The median corresponds to 

Qcrit ~ 2.5. 



B2 The midplane hydrostatic pressure of disk galaxies 

Under the assumptions of local isothermal stellar and gas layers, 
and a* > <t„ 



the midplane hydrosta tic pressure in di sks, P c: 



can be approximated to within 10% by dElmegreen|[l993l) 



— GE g 

2 s 



^Jgas H~~ ( " ) 



(B4) 



where E gas and E* are the total surface densities of gas and stars, 
respectively, and cr g and <r* give the vertical velocity dispersion of 
the gas and stars. We assume a constant gas velocity dispersion, 
a g — lOkms -1 (see §2.2.1). By assuming that E* 3> E gas , 
cr* = V 7rG h*E*, where h* is the stellar scale height. This ap- 
proximation could break down for very high redshift galaxies, 
whose disks are gas dominated. In such cases, we assume a floor 
of cr* ^ <j g . We estimate the stellar scaleheight assuming that it 
is proportional to the radial scalelength of the disk, as observed i n 
local spiral galaxies, with r e fj/ru ~ 7.3 ± 1.2 ( Krege fet al.l2002l) . 



B3 The KMT star formation law 

KMT theoretical model attempts to deal with three problems: (i) 
the determination of the fraction of gas in the molecular phase; (ii) 
the estimation of the characteristic properties (masses and surface 
densities) of GMCs, using a mixture of theoretical ideas and ob- 
served correlations; and (iii) the estimation of the rate at which 
molecular clouds convert themselves into stars. The latter rate 
is known observationally to be very small (~ 1% of mass per 
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free-fall time: iKrumholz & Tarj|2007h and is understood as a re- 
sult of the regulation of the SF activity by supersonic turbulence 
dKrumholz & McKedEjoll . KMT09 write the SFR per unit total 
gas mass as a function of two factors, the fraction of gas in the 
molecular phase, / mo j (see previous subsection), and the SFR per 
unit molecular mass, Vsf, as in Eq.(8] 

B3.1 The molecular-to-total gas ratio 

The fraction / mo i depends on metallicity and the gas surface den- 
sity of the atomic-molecular complex. The KMT09 expression for 
the molecular-to-total gas fraction predicts that in very metal poor 
environments, / mo i — > 0. If implemented at face value, this would 
prevent any SF from taking place in pristine gas at very high red- 
shift, and consequently there would be no enrichment to make pos- 
sible any subsequent SF. However, this behaviour results from ne- 
glecting gas-phase reactions for the formation of H2 molecules, 
which dominate over formation on dust grains when the metallicity 
is very low. Guided b y the results of calcula tions of the formation 
of the first stars (see lBromm & L arson 2004 for a review), we will 
therefore assume that a minimum = 10~ 4 applies at the onset 
of the SF activity in our implementation. 

B3.2 The star formation timescale 

KMT define the inverse of the timescale needed to consume the gas 
in a cloud into stars as ^sp = e ff/~tff, where tff is the freefall 
time of a GMC a nd e f f is the fraction of ga s converted to stars per 
freefall time (see Rrumholz & McKee 2005). In this model, e is a 
weak function of cloud properties, but the freefall time tff depends 
on the cloud mass M c i and surface density E c ;. KMT assume that 
the GMC mass corresponds to the critical Jeans mass for gravita- 
tional instability, and that GMCs are embed ded in a gaseous disk 
that is marginally stable by the lToomrd ( 119641) condition (Q = 1). 
KMT09 then use an observed correlation found for nearby galaxies 
between the global angular velocity of rotation around a galaxy and 
the gas surface density, which leads to M c i oc E gas . 



APPENDIX C: AN ILLUSTRATION OF THE IMPACT OF 
APPLYING DIFFERENT STAR FORMATION LAWS 

The semi-analytical model follows a range of physical processes, 
as set out in §2. The interplay between these makes it difficult to 
disentangle the impact of one ingredient in isolation. To gain an un- 
derstanding of the consequences of changing the SF law, we take 
the outputs of the original Bau05 and Bow06 models, and calcu- 
late the SFR using the SF prescriptions in Sj2] and compare with the 
SFR in the original model. To do this, we freeze the physical prop- 
erties of the galaxy used in the SFR calculation (i.e. galaxy size, 
stellar mass in the bulge and disk, cold gas mass, rate of accretion 
onto the disk of newly cooled gas, and metallicity of the cold gas) 
to make a consistent comparison, in which any difference will be 
due exclusively to the new SF law. Fig. lCll shows the ratio between 
the SFR calculated using a given SF law and that from the origi- 
nal recipe, SFRNcw/SFR or i g inai, as a function of stellar mass (left 
panel) and cold gas mass (right panel), at z — and z = 6, for the 
Bau05 and Bow06 models. 

SFRNcw/SFR or igi na i vs. stellar mass. At z = 0, the new SF 
laws result in larger SFR for low stellar mass central galaxies by 
around an order of magnitude. The exception is the KMT law in 
the case of the Bow06 model, which predicts a lower SFR than 



in the original model for low stellar mass galaxies, before rising 
and peaking around Mgtellar ~ 10 10 h -1 Mq. This arises from the 
abrupt fall in the radial Ssfr, predicted by the KMT law due to 
the large disk sizes of low mass galaxies in the Bow06 model (see 
Appendix D2). Moving to high stellar masses, in all cases the SFRs 
from the new laws are smaller than the original ones. In the Bau05 
model, similar trends are seen for satellite galaxies. However, the 
Bow06 model the predictions for satellites are quite different. This 
is because in the Bow06 model, satellite galaxies have, in general, 
more modest cold gas reservoirs compared with the Bau05 model. 
At z — 6, the predicted changes in the SFR are qualitatively similar 
to those seen at z = 0. 

SFRNcw/SFR or igi na i vs. cold gas mass. The largest differ- 
ences between the new and old SF laws are observed at interme- 
diate cold gas masses, 8 < log(M col d/h -1 Mq) < 10, while 
the differences become generally smaller when moving to the ex- 
tremes of high and low cold gas masses. This behaviour is more 
pronounced in the Bau05 model at all redshifts than in the Bow06 
model. The non-linear dependence of the SF laws on cold gas mass 
is directly linked with the non-linearity of the SF timescale with 
cold gas mass. Thus, we expect that changing the SF law also af- 
fects galaxies in different ways depending on their cold gas content. 



APPENDIX D: OTHER OBSERVED PROPERTIES OF 
GALAXIES 

We compare the predictions of the models with the new SF laws 
against selected observations. 

Dl The galaxy luminosity function 

The galaxy LFs in the 6j- and A"-bands at z — are shown in 
Fig. ID II for the Bau05 (top) and the Bow06 (bottom) models, and 
for variants with new SF laws. Perhaps surprisingly, nearly all of 
the new SF laws tested give a reasonably good fit to the observed 
LFs and follow closely the predictions of the original GALFORM 
models, even though we keep other parameter values fixed. Excep- 
tions arise for the KS and the KS. thresh laws in the Bow06 model, 
which produce too many bright galaxies in the blue band. However, 
the predicted A"-band LF is consistent with the observations even 
in these cases. This implies that the new laws result in too much SF 
activity at the present day in galaxies with M* 10 10 h _1 Mq for 
the KS SF laws in the Bow06 model. This suggests that the AGN 
feedback in these objects is not strong enough, and there is still gas 
cooling and hence SF activity at low-z. This lack of AGN feedback 
is a consequence of the reduced burst SF in these models at z ^ 4 
(see Fig. QJ, which reduces the mass of the BHs formed in mas- 
sive galaxies and therefore the strength of the AGN feedback (§2). 
Another indication of the absence of AGN feedback in the galaxies 
producing the excess in the fej-band, is that they correspond to late 
type galaxies with very low bulge-to-total stellar mass ratios (that 
host modest mass BHs). 

The changes in the LF are particularly mild in the Bow06 
model on applying the BR, BR.nonlin and KMT laws. This sim- 
ilarity is mainly caused by self regulation of the SF channels (i.e. 
starburst and quiescent SF), in which reduced quiescent SF activity 
at high redshift leads to galaxies characterised by massive disks and 
high gas fractions which experience more prominent SB activity in 
gas-rich mergers and disk instabilities. Thus the overall SF remains 
approximately constant. This is also reflected in the mild impact on 
the optical colours of galaxies, e.g. g — r, in which the blue cloud 
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Figure CI. Ratio between the new SFR and that in the original model, SFR^ CV ,/ S F R oligina i, plotted as a function of stellar mass (left-hand four panels) 
and of cold gas mass (right-hand four panels). The upper row shows the predictions at z = and the lower row shows z = 6. Within each set of four panels, 
the left column shows the Bau05 model and the right column shows the Bow06 model. Three representative SF laws are shown, the KS law (red), BR (blue) 
and KMT (green). Solid lines show the results for central galaxies only, while dashed lines show the results for satellite galaxies only. Only galaxies with 
A^cold > are included in this plot. Error bars show the 10 and 90 percentiles of the distribution, and for clarity are shown only for the KS law. 
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Figure Dl. The bj- (left) and if-band (right) galaxy luminosity functions 
at z = 0. The fiducial Bau05 (top) and Bow()6 (bottom) models are shown 
as black lines. Predictions are shown for the KS law (solid red), KS.thre sh 
law (dashed red; K98), BR law (solid blue; iBlitz & Rosolowskvl feOOol) ). 
BR.no nlin law (dashed blue; button et alj j20ld) ) and the K MT law (solid 
green; ItCrumholz et al. || 2009bl)). Observ ation al results from Norbe rg et alj 
|2002|) (6 T -band) and ICole et all feOOll) and lKochanek et alj <200ll) (K- 
band) are shown using grey symbols. 



in the colour-magnitude diagram appears slightly more pronounced 
when the new SF laws are applied. 

At high redshifts the similarity in the predicted LFs at z — 
is maintained. Fig. ID2I shows the rest-frame if -band LF at z — 1 
and z — 2 for the original models and variants. With the new SF 
laws the Bau05 model produces more galaxies at intermediate lu- 
minosities (i.e. — 16 < Mk — 51ogh < —20) at z — 2 compared 
to the original model. This arises because the bulk of the SF activ- 
ity in galaxies is shifted to higher redshifts with the new SF laws 
(Fig-UJ. Note that, even though the LF remains nearly unchanged, 
the contributions to the total luminosity from the bulge and disk 
components of galaxies are modified due to the change in the qui- 
escent and burst SF modes. The new SF laws tend to produce twice 
as many ellipticals as in the original models over the whole mass 
range. 

D2 Galaxy sizes 

In GALFORM the size of a galactic disk is determined by the con- 
servation of angular momentum of the gas cooling from the halo 
and the application of centrifugal equilibrium in the combined po- 
tential of the disk, bulge and host halo. Newly cooled gas modifies 
the angular momentum of the existing disk. Fig. lD3l shows the size- 
luminosity relation of late- (top) and early-type (bottom) galaxies 
in the two models (the Bau05 model on the left and Bow06 on 
the right), and the variants. Early-type galaxies are defined here 
as those with a dust extincted bulge-to- total luminosity ratio in 
the r-band exceeding 0.5 (B/T > 0.5; iGonzalez et al][2009h to 
approximately correspond to what is used for the SDSS samples. 
We also show two di fferent observational estimates of th e size- 
luminosity relations JShen et alj 120031 ; IDutton et all 1201 ll) . both 
based on SDSS data. 
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Figure D2. Rest-frame K -band galaxy luminosity function at redshifts z = 
1 and z = 2 as labelled, for the Bau()5 (left) and Bow06 (right) models and 
the different forms of SF laws tested: KS (solid red), KS. thresh (dashed red), 
BR (solid blue), B R.nonlin (dashed blue ) and KMT (solid g r een). Observa- 
tional resu lts from |Pozzetti et a l. 12003), Drory et al. 12004), Sarac co et alj 
l200el) and lCaputi et aiT l2006) are shown as grey symbols, identified by the 
key in the top-left panel. 



Considering first late-type galaxies, we find that the new SF 
laws produce only a small change in sizes for the Bau05 model, 
but a large change in the sizes of bright disk-dominated galaxies in 
the Bow06 model, improving the agreement with observations. At 
intermediate luminosities, both models agree somewhat better with 
the siz e-luminosity rela tion found by Duttor Tet al.1 J201 lh than with 
that o f lShenetal.l ( [2003h . The larger sizes reported bv lDutton et al.l 
are due to the 2-D fitting they perform to galaxy surface brightness 
profiles (i.e. including separate disk and bulge components, and in- 
clinations), in contrast with t he use of ci rcular apertures to estimate 
half-light radii in the case of lShen et all 

Bright spirals are larger in the Bow06 variants because the 
lower quiescent SFRs predicted by the new SF laws (§3.2) result 
in more massive disks compared to the original Bow06 model, 
which then causes more disk instabilities. (A lower SFR means 
less SNe feedback and consequently less cold gas ejected from 
the disk, hence giving a larger total mass - stars plus cold gas - 
for the disk). In the original Bow06 model, bright late-type galax- 
ies typically had quite large bulge-to-total luminosity ratios, close 
to the classification boundary B/T = 0.5, but the increased in- 
cidence of disk instabilities converts some of these to early-type 
galaxies. The late-type galaxies which remain have, on average, 
larger disk sizes than before (since smaller disks are more unstable 
for a given mass). In the Bau05 model, the differences are smaller 
since this model already produces longer SF timescales compared 
to the Bow06 model, due to the form and parameters used in the SF 
recipe (see §2.1). The failure of the Bow06 model to predict small 
enough sizes for faint disk-dominated galaxies might be amelio- 
rated with a more deta iled modelling of the radi al distributions of 
stars and gas (e.g. as in I Strin ger & Benson 2 0071) . 

For the early-type galaxies, on the other hand, the new SF 
laws result in large changes in sizes for bright spheroid-dominated 
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Figure D3. Galaxy half light radius as a function of r — band magnitude for 
late- (top) and early-type (bottom) galaxies. Early-types are defined to have 
a bulge-to-total luminosity ratio in the r-band exceeding 0.5. Different lines 
and colours show the medians in the mode l s teste d as indicat ed by the key . 
The observational results from Shen et al. 1 2003) and Dutt on et al] i201lh 
of late- and early-type galaxies (classified based on their Sersic index and 
colours, respectively), which are both based on the SDSS, are shown as 
grey symbols with errorbars (corresponding to the medians and 10%-90% 
of the distributions). The widths of the distributions around the median in 
the models are not shown for clarity, but are comparable to those of the 
observations (see Gonzalez et al. 20o3). 



galaxies in the Bau05 model, but only modest changes in the 
Bow06 model. The sizes of spheroids are driven by their forma- 
tion through mergers and disk instabilities. The flatness of the pre- 
dicted size-luminosity relations for spheroid-dominated galaxies 
in comparison with observations may result from the simplified 
calculation of the sizes of merger remnants, rather than implying 
that an additional me chanism is needed (see I Almeida etal. 2007; 
iGonzalez et alj2009l) . 



